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ABSTRACT 


An  attempt  Is  made  to  extend  ship-maneuvering  analysis  for  ap¬ 
plication  at  low  ship  speeds.  Captive-model  experiments  are  made  to 
Investigate  effects  of  non-equilibrium  propeller  speeds  on  rudder 
force  and  effective  thrust.  Other  experiments  are  made  to  explore 
first-order  effects  of  very  large  drift  angles  and  pure  yaw  rotation 
on  hull  hydrodynamic  reactions.  The  experimental  results ,  together 
with  previous  rotatlng-arm  data*  are  applied  in  motion  equations*  ac¬ 
counting  for  both  extreme  propeller  actions  and  the  occurrence  of 
large  drift  angles  and  turning  rates.  Also  included  are  propulsion 
and  rudder  time  lags,  and  the  Influence  of  uniform  water  current  on 
ship  trajectory.  The  equations  are  programned  for  computation*  and 
ship  response  to  a  simple  docking  maneuver  is  predicted.  The  compu¬ 
tation  is  repeated  several  times  to  examine  sensitivity  of  ship  mo¬ 
tions  and  trajectory  to  variations  of  operating  parameters*  namely* 
maximum  reverse  propeller  speed*  engine  response  time,  and  water  cur¬ 
rent  direction. 
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NOMENCLATURE 

The  basic  nomenclature  follows  that  of  SNAME  Technical  and  Research 
Bulletin  No.  1-5*  Principal  exceptions  are  the  choice  of  typical  dimen¬ 
sions  used  in  forming  dimensionless  quantities,  and  the  use  of  additional 
notation. 

Three  different  right-handed,  rectangular  coordinate  frames  are 

used: 

1.  (x.  y.  2)  are  body  axes  with  origin  at  the  center  of  gravity,  x  di¬ 
rected  toward  the  bow  along  the  longitudinal  centerline  hull  axis, 

y  directed  toward  starboard,  and  2  from  deck  to  keel.  The  x,  y 
plane  is  horizontal  with  the  ship  in  static  equilibrium. 

2.  (x  ,  y  ,  z  )  are  fixed  relative  to  the  earth,  with  the  x  and  y 

000  0  0 

axes  in  a  hori2onta1  plane  and  zp  directed  downward.  When  yaw 

angle  ♦  is  2ero,  the  x  and  y  axes  are  parallel  to  the  x  and 

0  0 

y  axes,  respectively. 

3.  (Xj,  y j ,  2j)  are  fixed  relative  to  the  body  of  water,  which  moves 
with  a  constant  velocity  relative  to  (x  ,  y  ,z  ).  This  velocity  is 

o  j  00 

called  the  water  current  and  is  denoted  IU  .  The  xT»  y.»  and  z, 
axes  are  always  parallel  to  the  xq,  y q,  and  z^  axes,  respectively. 

Unless  otherwise  noted,  all  dimensionless  hydrodynamic  (’hull)  deriva¬ 
tives  are  evaluated  at  zero  angular  velocity  and  acceleration,  zero  angles 
of  sideslip  and  rudder,  and  propeller  speed  set  for  model  self-propulsion 
point.  Definitions  of  symbols  used  are  listed  below: 

A  ,A1 , . .  • ,C  ,  etc.  Hydrodynamic  force  and  moment  coefficients 

'  used  in  motion  equations,  formed  from  products 

of  dimensionless  coefficients  and  dimensional 
quantities  p,  L,  H;  for  example,  Bj  =  p/2  LHY^ 

a  ,aa,..,b  ,bj  ,. .  ,c  ,  etc.  Dimensionless  coefficients  used  in  representa- 
*  “*  tions  of  propeller  and  rudder  forces 

CG  Center  of  gravity 

c^  A  propeller  thrust  coefficient,  useful  for  stalled  propeller 
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Propeller  diameter 

Bounding  values  of  l/Js»  distinguishing  different  regions  of 
propeller  operation}  used  for  selecting  relevant  representa¬ 
tion  during  computation 

Vector  sum  of  external  forces  acting  on  ship 

Ship  draft 

Ship  inertia  tensor 

Moment  of  inertia  of  ship  about,  z  axis 

Propeller  advance  ratio  in  behind-ship  condition 

Functions  used  when  obtaining  propeller  and  rudder  force  and 
moment  representat ions;  dimensionless 

Ship  length,  between  perpendiculars 

Mass  of  ship 

Hydrodynamic  moment  component  relative  to  z  axis;  yaw  moment 
Vector  sum  of  external  moments  acting  on  ship 
Propeller  speed  of  rotation 
Propeller  pitch 

Angular  velocity  component  of  ship  axes  in  z  direction 
Dimensionless  angular  velocity;  r*  =  rl/l) 

Limiting  value  of  r‘  for  crui si ng-type  maneuver 
Angular  acceleration  component  of  ship  axes  in  z  direction 
Time  constant  for  propeller-speed  response 
Time 

Apparent  slip  ratio  of  propeller 
Distance  along  trajectory 
Speed  of  CG  relative  to  fluid 
Speed  of  water  current  relative  to  earth 

Velocity  components  of  CG  relative  to  fluid,  measured  in 
directions  of  x  and  y  axes,  respectively 

Velocity  components  of  CG  relative  to  fluid,  measured  in 
directions  of  xq  and  yQ  axes,  respectively;  (Xj,y^)  and 

(xQ,yo)  coordinate  frames  are  always  parallel 

Velocity  components  of  CG  relative  to  earth,  measured  in 

directions  of  x  and  y  axes,  respectively 

o  o 

Velocity  components  of  water-current  relative  to  earth,  meas¬ 
ured  in  directions  of  xq,  yQ  axes,  respectively 

Acceleration  components  of  CG  in  directions  of  x  and  y 
axes,  respectively 


v  i  i  i 


x.v 


V'.»' 


yr,Nr 


*  V 

Hydrodynamic  force  components  relative  to  x  and  y  axes* 
respectively}  longitudinal  and  lateral  forces 

Typical  dimensionless  first  derivatives  of  a  hydrodynamic 
force  and  moment  with  respect  to  a  velocity  component 

Typical  dimensionless  first  derivatives  of  a  hydrodynamic 
force  and  moment  with  respect  to  an  angular  velocity  component 


Yl  ,NI 
v  v 


Typical  dimensionless  first  derivatives  of  a  hydrodynamic 
force  and  moment  with  respect  to  a  linear  acceleration  com¬ 
ponent 

Typical  dimensionless  first  derivatives  of  a  hydrodynamic 
force  and  moment  with  respect  to  an  angular  acceleration 
component 

Typical  dimensionless  third  derivatives  of  a  hydrodynamic 
force  and  moment  with  respect  to  an  angular  velocity  com¬ 
ponent  and  a  lateral  velocity  component 

X#Y»N(Hu11)  Sums  of  X>  Y ,  N  attributed  to  hull  motions*  respectively 


Yl  *NI 
r  r 


Y  '  «N 
rvv 


rvv 


X.Y.N 

(Prop.  +  Rudder) 


Sums  of  X.Y.N  attributed  to  propeller  and  rudder  actions 


P  Drift-angle;  sinp  =  -v/u 

P^  Limiting  drift  angle  for  crui si ng-type  maneuvers 

6  Rudder  deflection 

p  Hass  density  of  water 

i  Yaw  angle*  measured  about  z  axis*  relative  to  x  axis 

o  o 

tc  Direction  of  water-current  flow*  measured  about  z  axis 

relative  to  x  axis  ° 

o 

ui)  Angular  velocity  of  ship  relative  to  (xo*yo*zo)  axes 


Subscripts 

o  Referred  to  axes  fixed  with  respect  to  earth;  also  initial 

time  when  used  with  t 

1  Previous  steady  value 

*  Desired  or  commanded  value 


INTRODUCTION 

Upon  entering  pilot  waters  from  seaward,  the  difficulties  in  con* 
trolling  a  ship  increase  enormously  because  of  the  proximity  of  many 
solid  boundaries  such  as  harbour  Instal lat ions  and  other  ships.  The 
extremely  low  control  force  to  momentum  force  ratios  of  large  merchant 
ships  is  the  main  factor  which  makes  close  quarter  maneuvering  difficult. 

To  minimize  both  the  danger  and  damage  of  collision,  it  is  necessary  to 
drastically  reduce  ship  speed  in  this  type  of  operations. 

However,  the  reduction  of  ship  speed  introduces  other  problems 
which  complicate  ship  control.  Among  these  are  the  effects  of  wind  and 
current,  and  the  strong  dependence  of  rudder  foxes  on  propeller  speed 
(for  single  screw  ships).  All  these  factors  combine  to  cause  frequent 
saturation  of  propeller  and  rudder.  This  accounts  for  the  common  use 
of  tugs,  bow  thrusters,  or  anchors  during  close  quarter  maneuvers  of 
ships. 

The  work  done  under  the  present  task  order  was  divided  Into  two 

phases.  Ir.  the  first  phase,  the  effects  of  wind  on  low-speed  ship  con- 

2 ! 

trol  and  stability  were  treated  and  reported  separately  by  Eda.  The 
second  phase,  reported  herein,  deals  with  the  effects  of  steady  water  | 

currents  and  changing  propeller  speed.  In  neither  case  were  the  effects  , 

of  the  proximity  of  solid  boundaries  on  the  hydrodynamic  forces  treated. 

The  simplified  motion  equations  used  in  most  studies  of  ship  con¬ 
trol  provide  satisfactory  predictions  of  maneuvers  under  cruising  (con¬ 
stant  throttle)  conditions,  but  are  not  sufficiently  general  to  treat 
many  real  control  problems,  such  as  those  which  occur  at  low-ship  speeds. 

1  2 

Controllability  surveys  by  Norrbin,  Eda  and  Crane,  Williams  and 

34  5 

Noble,  Goodman,  and  Abkowitz  discuss  the  usual  forms  of  equations  and 
their  solutions.  None  Include  propeller-speed  changes  or  water-current 
effects.  Saunders,^  and  more  recently  Hawkins, ^  have  discussed  special 
low-speed  maneuvers  without  using  equations  of  motion.  Crane,  Uram  and 
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Chey  have  formulated  motion  equations  for  mooring  and  docking  maneuvers 
of  a  destroyer  and  a  submarine.  But  most  efforts  on  non-equilTbrlum 
proeel ler  effects  Have  concentrated  on  ship-stopping.  These  works  are 


summarized  In  Reference  9. 
II 


Papers  by  English,*®  Hawk  I ns, ^  Stuntz  and 


Taylor,' '  and  others  treat  bow-thrusters  and  other  maneuver  I ng-propuls Ion 
devices  (HPO) ,  but  these  do  not  include  HPD  forces  In  maneuvering  analy¬ 
sis. 


The  present  work  is  an  attempt  to  extend  the  generality  of  ship- 
maneuvering  analysis  to  include  engine -maneuver  and  steady  water-current 
effects,  and  to  introduce  the  possibilities  of  very  large  drift  angles 
and  dimensionless  turning  rates. 

Equations  of  motion  are  formulated  for  translation  and  rotation  of 
a  surface  ship  maneuvering  in  <;  reference  frame  whic!  is  fixed  with 
respect  to  a  moving  water  mass  (the  water  current).  The  independent  ship 
variables  comprise  propeller  speed  and  rudder  angle.  The  equations  are 
solved  with  respect  to  the  coordinate  frame  moving  with  the  water,  and 
the  solution  is  transformed  to  a  fixed  coordinate  frame.  Resultant  ship 
velocity  Is  integrated  to  obtain  the  position  and  orientation  of  the  ship 
at  any  time.  Water  current  is  a  parameter,  and  is  made  zero  when  dis¬ 
cussing  other  effects. 

Coefficients  of  the  equations  are  obtained  f rom  exper iment ,  cal¬ 
culation,  and  existing  data.  The  combined  effects  on  the  hydrodynamic 
forces  and  moments  of  ship  speed,  propeller  speed,  and  rudder  angle 
(at  maneuvering  speeds)  are  investigated  by  captive-model  tests  in  a 
straight  towing  tank.  Indications  of  hydrodynamic  reactions  to  very 
large  drift  angles  and  to  pure  yaw  rotation  are  obtained  experimentally. 
Numerical  integrations  are  performed  by  digital  computer. 

A  few  examples  of  computed  maneuvers  are  given  to  illustrate  uses 
which  can  be  made  of  more  general  analyses  of  ship  maneuvering.  Using 
the  computation  procedure  developed,  the  sensitivity  of  ship  response 
to  particular  ship  parameters  is  studied  for  a  simple  ship-docking 
maneuver. 

It  is  noted  that  some  of  the  procedures  used  in  formulating  the 
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equal l  ons  of  motion  are  merely  expedient  measures  taken  to  permit  qualita¬ 
tive  results  to  be  obtained.  More  fundamental  treatment  of  these  close- 
quarter.  low-speed  control  problems  Is  left  for  future  studies. 


This  project  was  sponsored  by  the  Bureau  of  Ships  General  Hydro¬ 
mechanics  Research  Program  under  Contract  Nonr  263(63)  and  technically 
administered  by  the  Oavid  Taylor  Model  Basin. 


EQUATIONS  OF  MOTION  -  GENERAL 


General  Form  of  Motion  Equations 


Differential  equations  describing  the  motion  of  a  ship  are  written 


with  respect  to  a  reference  frame  rotating  with  the  ship's  axes. the  origin 
of  which  is  chosen  coincident  with  ship's  center  of  mass  (see  Figure  ’)• 

F  =  m(^  +  x  V) 

< 1 ) 

7?  =  I  +  uj  x  r  533* 


where 


F  sum  of  all  externa!  forces 

N  sum  of  all  externa!  moments 

V  velocity  of  ship  center-of-mass 
to  angular  velocity  of  ship  (ship's  axes) 
I  moment  of  inertia  tensor 


All  vectors  can  be  expressed  in  rotating  (ship  axes)  coordinates. 

For  a  surface  ship  in  calm  water,  hese  vector  equations  are  reduced  to 
three  Cartesian  equations  by  applying  the  following  simplifying  assump¬ 
tions: 

1.  Ship  is  constrained  to  move  in  its  lateral  plane 
(i.e.»  the  x,  y-pl ane) * whi ch  is  horizontal. 

2.  Center-of-mass  is  in  the  principal  plane  of  symmetry  of  the  hull. 

3.  Mass  and  center  of  mass  are  constant. 

4.  Products  of  inertia  are  zero. 


Then , 

m(u  -  rv)  =  X 

m(v  +  ru)  =  Y  (2) 

I  r  =  N 

z 

where  X  and  Y  represent  summations  of  unbalanced  external  forces  act¬ 
ing  in  the  x  and  y  directions,  and  N  the  summation  of  unbalanced 
moments  about  the  z-axis. 
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External  forces  arise  from  hydrodynamic  reactions  to  hull  motion, 
propeller  and  rudder  actions,  aerodynamic  reactions  on  hull  structure, 
and  all  other  forces  which  may  be  applied  by  waves,  bow-thruster,  anchor 
chain,  tugboat,  mooring  lines,  etc.  In  this  analysis,  only  hull  hydro- 
dynamic  and  propeller-rudder  forces  are  included.  In  succeeding  sections, 
their  contributions  to  X,  Y, and  N  are  considered. 

Reference  Frames 

It  is  clearly  desirable  to  describe  low-speed  ship  maneuvers  relative 
to  an  inertial  reference  frame  fixed  in  the  earth.  However,  the  motion 
equations  and  solution  are  greatly  simplified  if  referred  to  a  coordinate 

V: 

frame  fixed  in  the  ship,  with  subsequent  transformation  of  solutions  to 
the  desired  reference  frame  (  Figure  ll 

With  a  uniform  water  current,  an  additional  (intermediate)  reference 
frame  is  useful.  The  additional  frame  maintains  the  same  orientation  as 
the  fixed  frame,  but  moves  with  the  steady  velocity  of  the  uniform  water 
current.  In  computation,  the  first  transformation  of  the  solution  is 
from  ship  axes  to  intei  .nedi ate  axes.  The  second  transformation  simply 
adds  the  water-current  velocity  components,  to  obtain  motions  relative  to 
the  final  desired  reference  frame. 

It  is  emphasized  that  the  uniform  water  current  need  not  enter  the 
differential  equations  of  motion,  because  a  steady  velocity  will  intro¬ 
duce  no  acceleration  in  the  motion  equations.  The  current  is  entirely 
accounted  for  in  the  final  integration  of  velocities.  Steps  involved  in 
the  velocity  transformations  are  shown  below 

(1)  Ship  velocity  components  (expressed  in  ship  axes)  are  transformed  to 
intermediate  axes  by  a  rotation  of  coordinates. 

Uj  =  u  cosilr  -  v  sin^ 

(3) 

Vj  =  u  s i n\lf  +  v  cos^ 

Orientation  of  ship  axes  with  respect  to  ship  geometry  varies  according 
to  ship  loading  conditions  because,  for  convenience,  x  and  y  ship-axes 
are  taken  horizontal  (hence,  z-axis  vertical). 
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and 


rd  t  +  t(o) 


where 

u  and  v  are  coordinates  of  velocity  along  ship  axes*  and 
Uj  and  Vj  are  components  of  ship  velocity  along  intermediate  axes 


(2)  Components  of  water-current  velocity  ari  then  added  to  obtain  ship 
motion  relative  to  fixed  inertial  axes. 


u 

o 


+ 


v 

o 


+  V 


or 

u  =  u  cos^r  -  v  sin^r  +  U  cos^r 

o 

(4) 

v  =  u  sin^r  +  v  cos^  +  U  sir.^r 

o 

keeping  in  mind  that  the  intermediate  axes  (I)  and  fixed  inertial  axes 
(o)  are  paral lei . 


Time  integration  of  the  velocities  uq  and  vq  yields  instantaneous 
ship  position  ( xq  ,  yQ)  vith  respect  to  any  chosen  origin. 
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HULL  CONTRIBUTIONS  TO  EXTERNAL  FORCE  AND  MOMENT 


Hydrodynamic  forces  exerted  on  a  ship's  hull  are  attributed  to  vis¬ 
cous  and  inertial  fluid  effects  which  are  expressed  as  functions  of  hull 
motions  and  accelerations.  Since  a  wide  range  of  maneuvering  motions  may 
occur  at  low  ship  speeds,  hydrodynamic  representations  are  separated  ac¬ 
cording  to  application  to  crui sing-type  or  docking-type  maneuvers.  The 
former  case  is  considered  first. 

A  convenient  form  for  expressing  hull  hydrodynamic  forces  is  avail- 

1  2 

able  in  the  Taylor  series  expansion.  In  the  present  analysis,  the 
general  expansion  is  specialized  to  the  surge,  sway,  and  yaw  motions  of 
a  ship. 

The  following  simplifications,  extra  to  those  noted  in  the  previous 
section  ,  are  made: 

1.  The  hydrodynamic  derivative  coefficients  X^  ,  X^J.  ,  Y£  and 

N^  are  assumed  to  be  negligible.  For  convenience  a  term  pro¬ 
portional  to  uU  is  used  to  account  for  X'  and  drift-angle 
effects.  ° 

2.  Hydrodynamic  terms  of  order  four  ard  higher  are  assumed  to  be 
negl igible. 

3.  Forces  due  to  hull  motions  and  accelerations  are  represented  to 
be  independent  of  propeller  and  rudder  actions. 

4.  Ship  floats  in  deep  ( O/H  >  3*5)»  calm,  unrestricted  water. 

Simplification  (3)  is  necessary  due  to  the  absence  of  captive-model  ex¬ 
perimental  data  for  a  ship  on  curved  path  with  reversed  propeller  rota¬ 
tion,  or  in  any  propeller  condition  far  removed  from  model  or  ship  self- 
propulsion  point.  An  indication  of  the  effects  of  non-equilibrium 
propeller  speed  on  hull  hydrodynamic  data  was  obtained  by  inspection  of 
rotating-arm  tests  results  of  a  model  of  T.S.  EMPIRE  STATE  IV. The 
measured  effects  cannot  be  explained  simply  on  the  basis  of  propeller- 
side-force  changes  at  different  propeller  loadings.  Although  the  differ¬ 
ences  are  large  enough  to  warrant  further  investigation,  they  are  assumed 
to  be  of  second-order  importance  for  present  purposes. 
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With  the  preceding  simplifications,  the  nonlinear  hull  for ce-and-mo- 
ment  expressions,  retaining  third-order  terms  in  the  velocities,  may  be 
written  as 

X,  ..  =  A  u  +  A.uU  *■  A «rv 
hull  o  1  2 

Yh"i1  =  V  +  Bluv  +  B2ur  +  B3rs^/U  +  8 4^ v/u  +  B^rVu  +  B^vVu  ►  (5) 

N.  , ,  =  C  r  +  C.uv  +  C-ur  +  C,rv®/u  +  C.  rav/u  +  C_r3/u  +C,va/u 
hull  o  1  2  3  4  '  5  6 


where  It  has  been  assumed  that  |  v/u  |  <  0.26  and  r'  <  1.0  ,  and 


A  =  ■?  LaHXl  B  =  -x  La  HY 1 
o  2  u  o  2  v 


A,  =  f  LHX' 

A2  =  f  L=HXri 


B.  =  £  LHY ' 

1  2  v 

B2  =  |  Ls  HY^ 


Co  =  \  L*HN. 

c,  =  f  l»hn; 

C2  =  |  L3  HN^ 


B,  =  £  LaH(Y  '  +  Y 1  )  C  =  £  L3H(N  '  +  N') 
3  4  '  rvv  r  34  rvv  r 


B.  =  77  l3HY  • 
4  4  rrv 


6S  =r2l,mrrr 


C4  *  5  L‘H  "rri 
C5  =  f2  L’H  N rrr 


B,  =  r>  LH(Y  '+3Y')  C,  =  p,  LEH(N  '  +  3N •  ) 
6  12  vvv  v  6  12  vvv  v 


All  dimensionless  coefficients  except  XI,  X  '  ,YI  ,  and  N1  are  evaluated 

r  u  vr  v  r 

experimentally,  as  described  in  Reference  13* 

The  derivatives  X^l  ,  Y^  ,  and  NI  relate  hydrodynamic  forces  to  body 

accelerations,  and  X  '  is  taken  equal  in  magnitude  to  Yl  .  The  Yl 

vr  v  14  v 

and  Nl  coefficients  are  computed  following  the  method  of  Lewis  and 
f  1 5 

Prohaska.  For  example, 

.  fXs 

7  ^C  ’K  k  (x^dx 

-Y./m  =  - - - - - - - - 

V  fix 

I  S  S(x)dx 
xb 


In  this  strip-wise  integration,  k^  is  the  two-dimensional  lateral  added- 
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mass  coefficient.  This  is  determined  for  ship-like  sections#  by  using 

Prohaska’s  results.  The  k  is  a  three-dimensional  correction  factor 

c 

determined  by  comparing  the  exact  added-mass  computation  of  the  prolate 
spheroid  given  by  Lamb  with  that  obtained  by  the  stripwise  computation 
for  the  same  spheroid.  Factor  S(x)  is  the  local  sectional  area.  The 
rotary  acceleration  derivative  N!  is  obtained  in  a  similar  manner,  but 
includes  an  x"3  factor  in  the  integrand  of  the  numerator.  The  longitudi¬ 
nal  accele’ation  derivative  X^l  is  taken  equal  to  that  of  an  equivalent 
prolate  spheroid. 


Conditions  at  Very  Low  Speeds 

Experience  shows  that  at  moderate  ship  speeds,  drift  angle,  p  •  does 
not  normally  exceed  approximately  ±1 0  degrees,  and  dimensionless  turning 
rate,  r1  ,  is  normally  limited  to  about  0.7  (corresponding  to  turning- 
diameter/ship-length  ft*  3).  However,  a  ship  maneuvering  at  very  low  speeds 
may  develop  much  larger  drift  angles  and  space  turning  rat3s. 

Consider  the  definition  of  P  ,  as  forward  speed  u  becomes  small 
with  respect  to  lateral  speed,  v  : 

lim  p  =  I  im  ["-tan  ^  (^)J  =  *  w/2 
u— o  u-o 

In  similar  fashion,  space  turning  rate  grows  Targe  as  ship  speed  vanishes; 

i  *6 •  i 

lim  r1  =  lim 
U-o  U-o 

Such  situations  arise  with  the  action  of  external  forces  which  are  not 
wholly  dependent  on  ship  motions. 

Two  broad  categories  of  maneuvers  are  defined: 

1.  Crui si ng-type  maneuvers:  Both  of  the  below  criteria  are  satis¬ 

fied. 

2.  Oocking-type  maneuvers:  One  of  the  below  criteria  is  un¬ 

satisfied,  or  both  are  unsatisfied. 

I  P!  <  PL  !  r*  !  <  r- 


9 


R - 1 1 69 


These  definitions  are  useful  in  defining  maneuvers  for  which  con¬ 
ventional  hydrodynamic  expansions  are  valid.  These  are  called  cruising- 
type  maneuvers#  and  it  is  for  these  that  most  experimental  data  are 
obtained.  Ranges  of  variables  may  be  extended  to  cover  docking-type 
maneuvers,  of  course  (given  sufficient  date),  but  it  is  possible  that  an 
expansion  other  than  Taylor's  series  may  be  more  convenient  for  treating 
docki ng-type  maneuver s  (double  Fourier  series,  for  example).  If  a  Taylor 
series  is  used,  an  expansion  point  at  (u=v=r =u=v=r =0)  might  prove  more 

convenient  than  at  (u=u  ,  v=r =u=v=r=0) . 

o 

In  this  report  an  example  of  a  simple  docking  maneuver  is  computed. 
The  maneuver  is  terminated  as  forward  motion  stops  after  sustained  re¬ 
versed  propeller  action.  The  computation  will  show  a  tendency  for  the 
ship  to  develop  large  values  of  p  and  r'  at  the  end  of  the  maneuver. 

To  allow  for  this,  terms  exhibiting  first-order  qualities  of  the  true 
(steady-state)  hydrodynamic  reactions  are  used  in  this  region.  Coefficients 
of  the  simplified  terms  are  relatively  easy  to  evaluate,  since  they  repre¬ 
sent  uncoupled  drift-angle  and  yaw-rate  effects.  The  acceleration  de¬ 
pendent  terms  appear  as  before.  Simplifications  3  and  4  still  apply. 

hull  o  1  2 


=  B  v  +  B_vU  +  Bfl  v  |  v|  +  BQur 


Nu  . .  =  C  r  +  C,vu  +  CQ  — —  I  vu  |  +  CQr  I  r| 
hull  o  7  8  us  1  '  9 

where  it  is  assumed  that  u  >  0  . 


The  relative  importance  of  steady-state  hydrodynamic  forces  obviously 
declines  (relative  to  inertia  effects)  as  ship  speed  decreases.  However, 
it  cannot  simply  be  stated  that  steady-state  terms  are  negligible  below 
any  particular  speed  (except  under  very  special  conditions).  For  example, 
consider  ship  response  to  a  pure  applied  yaw  moment.  After  initial  yaw 
acce I erat ion ,  the  resisting  hydrodynamic  moment  is  almost  entirely  a 
steady-state  damping  effect.  If  steady-state  terms  were  arbitrarily  ex¬ 
cluded,  acceleration  without  limit  would  be  computed. 
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PROPELLER  AND  RUDDER  CONTRIBUTIONS  TO  EXTERNAL  FORCE  AND  MOMENT 


An  important  departure  from  usual  maneuvering  analyses,  at  low  ship 
speeds,  is  the  loss  of  proportionality  between  control  forces  and  the  hull 
forces  which  are  associated  with  motion  through  the  water.  The  control 
forces  of  propeller  and  r  u  d  de  r  can  be  di  sproportionately  large  or 
small,  according  to  the  action  of  the  propeller.  External  forces  and  mo¬ 
ments  attributed  to  propeller  and  rudder  actions  depend  on  several  factors 
including 


1.  Propeller,  hull,  and  rudder  configurations 

2.  Propeller  speed  of  rotation  and  angular  acceleration 
3>  Ship  linear  and  angular  velocity 

4.  Rudder  angular  deflection  and  deflection  rate 

5*  State  of  cavitation  or  ventilation  of  propeller  and  rudder 

6.  Ship  hydrodynamic  environment,  i.e.,  water  density,  depth, and 
lateral  boundaries,  and  waves 


The 
simpl i fy 


following  assumptions,  together  with  those  in  previous  sections, 
representations  of  propeller  and  rudder  forces: 


1.  Ship  speed  is  low;  hence  lateral  ship  velocity  is  small. 

2.  Yaw  rotation  is  slow. 

3.  Based  on  ( 1 )  and  (2),  the  lateral  component  of  inflow  to  pro¬ 
peller  and  rudder  is  small. 

4.  Rudder  deflection  rate  is  slow,  conforming  to  normal  ship  prac¬ 
tice  (  2-l/2°/sec).  Also,  frequency  of  rudder  oscillation  is 
low;  hence  errors  due  to  quasi-steady  representation  are  small. 

5.  Accelerations  of  propeller  rotation,  rudder  deflection,  and  hull 
motion  which  might  affect  propeller  and  rudder  forces  have  negli¬ 
gibly  small  effects  on  integrated  hull  motions. 

6.  Propeller  and  rudder  are  not  cavitating. 


With  the  above  assumptions,  propeller  and  rudder  forces  are  repre¬ 
sented  as  functions  of  propeller  speed,  n  .ship's  longitudinal  speed,  u  , 
and  rudder  angular  deflection,  6  . 


1 1 


X  Force  Due  To  Propeller  and  Rudder 

The  dimensionless  thrust  coefficient,  k  ^  =  t^ir  JSt  ,  is  commonly  used 

t  pd4  na  7 

to  represent  propeller  test  results  for  various  conditions  of  ship  speed 
and  propeller  rotation.  This  notation  is  useful  as  an  aid  in  understand¬ 
ing  and  representing  propeller  force  effects.  In  this  study,  self-pro¬ 
pelled  captive-model  tests  have  been  used  to  measure  a  quantity  which 
will  be  denoted  X(prop.  rudder).  This  quantity  is  defined  as  total 
measured  axial  force,  X(u,n,6)«minus  the  resistance  of  the  model  hull 
without  propeller,  X(u,— , — ).  From  this,  a  dimensionless  "effective 
thrust"  coefficient  is  defined. 

k  _  X(u  ,n ,  6 )  -  X(u  ^ 

x  p  d4  n° 

As  with  open-water  propeller  thrust  characteristics,  it  is  convenient  to 

graph  k^  as  a  function  of  propeller  advance  ratio  (here  using  ship's 

speed  of  advance),  i.e..  J  =  •  For  clarity,  rudder  effects  are  o- 

s  nd  ' 

mitted  from  the  next  discussion. 

Sketch  I  shows  the  relation  of  k^  to  Js  in  regions  of  propeller 
operation  which  are  important  for  ship-stopping. 


The  path  of  kx  during  rapid  propeller-reversal  is  traced: 

1.  Approach  to  equilibrium  condi t  ion  (thrust  =  resistance) 

2.  Propeller  speed  reduced;  blade  sections  approach  zero  angle  of 
attack  (thrust  zero) 

3‘  Propeller  speed  further  reduced;  blade  sections  at  negative 
angles  of  attack  (thrust  negative) 

4.  Propeller  speed  further  reduced;  flow  about  blades  completely 
separates  (propeller  dragging) 

5*  Propeller  stopped  and  dragging*  ...  6.  Propeller  turning 

slowly  astern  but  dragging 

7*  Astern  propeller  speed  increased  (blades  recover  from  stall  and 
begin  thrusting  astern) 

8.  Propeller  speed  approaches  a  constant  value;  ship  speed  gradually 
reduces;  astern  thrust  determined  by  propeller  characteristics 


The  representation  of  Sketch  I  is  not  useful  in  the  interval  between 
4  and  7*  Mathematically,  as  n  approaches  zero.  tends  to  infinity. 
Furthermore,  beyond  the  blade  stall -point  at  large  negative  angle*  of 
attack,  k  cannot  be  represented  simpiy  in  terms  of  J  .  For  cargo 
ships  and  tankers  the  negative  blade-angle  stall  will  occur  at  a 
value  between  2  and  4.  Another  dimensionless  coefficient  is  then  con¬ 
venient  to  express  "effective  propeller  thrust."  This  is  the  coefficient 
X  1 6 

C*  =  p  \jP  '  USeC*  ky  B i nde  1  and  others.  The  singularity  in  Js  =  ^ 
is  avoided  by  using  the  reciprocal,  ^  . 


SKETCH  II 
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The  characteristic  behavior  of  c^  in  a  ship-stopping  maneuver  is 
shown  in  Sketch  II  (numbers  correspond  to  those  of  Sketch  I).  Because  pro¬ 
peller  speed  will  change  much  more  quickly  than  ship  speed*  the  behavior  of 
X(prop.)  is  similar  to  that  of  c^  in  Sketch  II. 

Because  of  slow  propeller  spteds  in  the  region  from  4  to  6,  experi¬ 
mental  evaluation  of  c^  is  subject  to  serious  scale  effect  ( especial 1y  for 
propeller  operating  behind  the  ship);butif  c^  passes  from  4  to  6  in  a  few 
seconds,  the  accuracy  of  the  representation  in  this  region  is  not  crucial. 

An  estimated  propeller  drag  coefficient  is  therefore  used  (See  Appendix  a)* 

When  unstalled  propeller  operation  resumes  (at  some  reversed  propeller 
speed),  the  k^  vs.  relation  is  again  effective  and  simple  to  use.  The 
required  segments  of  the  k^  ,  curve  c'*e  easily  described  by  simple  poly¬ 
nomials  in  Js  .  For  example,  in  the  rerion  from  1  to  4  (of  Sketch  I),  an 
adequate  curve-fi t  is  obtained  using  k^  =  aQ  +  aj  J&  +  ^  .  And  since 

X(prop)  =  kx  pd4  rp 

substitution  of  the  polynomial  in  Js  for  k^  *  and  for  *  yields 

X(prop)  =  (pd4aQ)na  +  (ptPa^un  +(pd2aa  )u3  (8) 

where  a^,  ax  ,  a3  are  determined  for  a  particular  propeller  and  hull. 

In  similar  fashion,  the  k^  curve  segment  passing  from  point  7  through 
point  8  may  be  curve  fitted  to  .  For  digital  computation  of  maneuvers, 
the  correct  expression  for  k^  can  be  automat i cal  1 y  selected  according  to 
the  i nstantaneous  value  of  ,  using  rules  based  on  Sketch  I. 

Rudder  effect  on  X  force  is  treated  next.  The  drag  of  a  rudder 
mounted  aft  of  an  ahead-turning  propeller  is  a  function  of  rudder  angle, 
ship  speed, and  propeller  speed.  If  k^  is  approximately  parabolic  with 
rudder  angle.  6  ,  a  simple  extension  of  the  propeller-thrust  representation 
to  rudder  drag  may  suffice. 

X(rudder)  ;  ^  pd*n®  6s  (9) 

Here  we  express  5  =  +  aA  Jg  ,  leading  to 
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X( rudder)  =  £(pd4  a,  )na  +  ( pd3  a4  )un  J  6a  (10) 

When  propeller  is  dragging  or  reversed*  the  rudder  drag  Is  small  and 
erratic.  No  attempt  is  made  to  represent  it. 

Y  Force  Due  to  Propeller  and  Rudder 

Representations  of  Y  force  and  N  moment  due  to  propeller  and 
rudder  closely  follow  those  for  X  force.  The  following  will  pertain  to 
both  Y  and  N  . 

Because  the  flow  past  the  rudder  (of  a  single-screw  ship)  is  de¬ 
termined  by  both  ship  speed  and  propeller  speed*  the  concept  of  k  vs.  J 
is  useful  for  developing  a  polynomial  representation  for  side  force*  but 
now  only  in  the  propeller  operating  region  from  1  to  4  of  Sketch  I.  If 
in  this  region  Y  is  a  linear  function  of  rudder  angle*  the  symmetrical 
part  of  Y(prop.  +  rudder)  with  respect  to  6  may  be  represented  by 

£(pd4  bo)na  +  (pd1  bj  )un  +  (pdaba)uaJs 

A  small  unsymmetrical  side  force  is  also  observed  for  single-screw  ships. 
This  is  called  the  Hovgaard  effect^  and  is  caused  by  asymmetry  of  propeller 
rotation.  Several  authors  attempt  to  explain  this  physically*  but  it  is 
represented  here  as  a  simple  function  of  u  and  n  ,  based  on  results  of 
straight  course  experiments  in  the  next  section.  Summing  the  symmetrical 
and  unsymmetrical  parts*  the  total  expression  is 

Y(prop.  +  rudder  )=£(pd  4bQ)na  +( pd1!^  )un+(  pdaba)ua  J  6+(pd4b3)  na+(  pdab4  )ua 

(ID 

In  the  stalled  region,  for  advance  ratios  to  the  right  of  point  4, 
flow  aft  of  the  propeller  is  much  disturbed.  Rudder  forces  are  then 
sharply  reduced  and  difficult  to  evaluate.  In  this  condition,  it  is  es¬ 
timated  that  rudder  effectiveness  is  reduced  to  less  than  l/3  of  that  in 
normal  ahead  operation.  This  is  based  on  model  tests  of  a  related  single¬ 
screw  hull  with  propeller  removed.^ 

As  propeller  rotation  is  reversed  the  propeller  race  is  directed 
forward.  Local  flow  past  the  rudder  is  destroyed  and  steering  control  is 
lost.  Side-force  bias  due  to  propeller  rotation  gains  importance*  and  it 
aione  influences  the  directional  behavior  of  the  ship  (in  the  absence  of 
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other  external  forces). 

N  Moment  Oue  to  Propeller  and  Rudder 

The  representation  of  yaw  moment  produced  by  propeller  and  rudder, 
N(prop. +  rudder  ) ,  closely  follows  that  for  side  force.  An  additional 
length  dimension  appears  in  the  exponent  of  d  .  The  number  of  terms  in 
will  depend  upon  che  shape  of  the  function  k^,  and  the  precision  re- 
qui red. 

A  summary  of  terms  used  in  the  present  study  for  N(prop.  + rudder ) 
appears  below: 

Operating  Region 

of  Propel  1 er  Expression  for  Rudder-PropeH er  Yaw  Moment 


to  4 

N,  =(pdsc  n®  +  pd  4c,  un  ♦  pd^c-  us)6  +  pdB  c,  n  ® 
c  ,n  '  o  1  * 

(12) 

to  6 

Nk  =  pd3  c.ua  6 
o  «n  * 

03) 

to  8 

N6.n= 

(14) 

Although  insufficient  data  are  available  to  analyze  ship  dynamic  stability 
in  this  condition,  the  reversed  propeller  may  destabilize  the  ship  by  its 
effect  on  the  flow  about  the  hull.  While  this  could  cause  yaw  divergence 
in  either  direction,  at  low  speeds  the  unbalanced  side  force  produced  by 
the  reversed  propeller  usually  causes  yaw  angle  to  de/elop  in  the  posi¬ 
tive  ( cl ockwi se ) sense.  This  tendency  is  shown  in  loter  computations.  The 
yaw  effect  is  used  to  advantage  by  ship-handlers  and  explains  their  prefer¬ 
ence  for  port-si de-to-pier  landings. 
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TIME  LAGS  OF  PROPELLER  ANO  RUOOER  RESPONSE 


Engine  Orders 


Commands  from  the  bridge  to  engine-control  indicate  desired  magnitude 
and  sense  of  propeller  rotation.  A  typical  engine-order  (bell)  table  for 
a  merchant  ship  is  shown: 


Engine  Order 

Full  Ahead 
Half  Ahead 
Slow  Ahead 
Dead  alow  Ahead 
Stop 

Slow  Astern 
Half  Astern 
Full  Astern 


Propeller  Speed  (rpm) 

60 

40 

20 

10 

0 

-15 

-30 

-45 


The  system  used  aboard  naval  ships  is  slightly  different.  In  the 
naval  system,  the  engine  (propeller)  is  specified,  followed  by  the  desired 
sense  and  speed  of  rotation  (in  thirds).  For  this  study  a  merchant  ship 
is  assumed. 


Engine  response  to  orders  requiring  large  amounts  of  power  is  limited 
by  the  steam  pressure  available.  Minimum  allowed  levels  are  established  to 
protect  the  boiler  against  excess  steam  drain  through  the  astern  turbine. 
These  limits  guide  throttle  opening  until  desired  propeller  speed  is  ap¬ 
proached. 


Bridge  control  is  not  treated  here,  but  may  be  introduced  immediately, 
since  there  is  no  restriction  to  discrete  propeller  speeds  in  this  analy¬ 
sis. 
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Propeller  Time  Lag 


Rate  of  change  of  propeller  speed  will  depend  upon 


1-  Instantaneous  propeller  speed,  n 
2.  Instantaneous  longitudinal  speed,  u 

3»  Previous  steady  propeller  speed,  n.  ,  and  nozzle  combination 
4.  Propeller-speed  ordered,  n 


A  complete  analysis  of  transient  propeller  speed  (including  pro¬ 
pulsion-machinery  characteristics  and  human-response  factors)are  beyond  the 
scope  of  the  present  study.  A  simplified  function  is  used,  which  provides 

a  good  approximation  for  the  purposes  of  this  work.  It  is  based  on  full- 

19 

scale  ship  data  for  the  "crash-back"  maneuver.  Propeller  response  is 
described  by  the  first-order  differential  equation  -g^  +  Cjn+Cg  =  0. 
Applying  boundary  conditions  n  =  nQ  at  t  =  tj  and  n  -*  n*  as  t  -*  <»  » 
instantaneous  propeller  speed  is  given  by 


(n*  - 


nf)(l  -e’At/T) 


At  =  t-t 


(15) 


where  T  =  time  constant  for  particular  ship  and  maneuver.  This  may  be 
estimated  by  fitting  Equation  (15)  to  given  response  data. 


Rudder  Time  Lag 


The  steering  machinery  of  a  large  ship  will  provide  approximately 
constant  rudder  rotation  For  rudder  changes  of  more  than  a  few  degrees. 
According  to  U.S.  Coast  Guard  Marine  Engineering  Regulations  and  American 
Bureau  of  Shipping  Rules,  average  rudder  rotation  shall  not  be  less  than 
2-1/3  degrees  per  second.  In  the  S.S.  GOPHER  MARINER  trials,  the  rudder 
was  shifted  from  35°  left  to  35°  right  in  22  seconds,  i.e.,  at  an  average 
rate  of  3«2  degrees  per  second. 


by 


For  constant  rudder  angular  rate,  rudder  deflection  is  represented 

•  '* 

6  =  6.  +  fi(t-t.),  and  6  >  6  according  to  sense  of  rotation,  where 


\ 
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6^  =  initial  deflection  angle 

tj  =  time  at  execution 
• 

6  =  rudder  angular  rate 

*iV 

6  =  ordered  rudder  angle 

•  Vf 

The  sign  of  5  is  determined  by  the  sign  of  (6  —  6.  )  .  Rudder-angle 
mands  are  specified  in  degrees  right  or  left  of  amidship  (-or  +  » 
spect i vel y  ) .  but  radian  measure  is  us^d  in  equations. 


corn- 
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SUMMARY  Or  MOTION  EQUATIONS 


The  equations  of  motion  are  summarized  below  In  a  form  for  solution! 


(m  -  A  )u  = 
o 


AjUU  ♦  (m  ♦  A^)rv 


X(Hu1 1  ) 


N  EQUATION 


Ajna+  A^un  *  A^u 3 ♦  (A ^n®1  ♦  A^un)8a  . 


*/Js  >  g, 


AguU  ♦  A^u9  K* 


A|0n*  *  A» |Un 


'  91  >  ,//js  >  g2 


.  g2  »  l/Js  >  g3 


X(P'op  ♦  Rudder) 


*1  2°*  *  A|  3un  *  A|4u:>  *A,5uVn 


9j  >  1/J„ 


(m  -  Bq)v  = 


B  jUV  ♦  ( B^-m^ur  ♦  Bjrv3/u  ♦  B^r’v/u  «  B^r-'/u  ♦  B^v^/u  l  P  <  p^ ,  r 1  s  r£,and  u  i  u^ 


Y( Hu  1 1  ) 


BjVU  ♦  Bgv|v|  ♦  (Bg-m)ur 


I  p  >  p^or  r‘>  r^  or  u  <  u^ 


B|0n3  +  +(B)2n3  *  BHun  "  9iOjB)6  '  l/Jf  >  g. 


M  »|5U*# 


»  g,>  »/Js>  g2 


,  g2.-  »/j# 


Y(Prop  ♦  Rudder ) 


(i,-co)M 
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,/Js  >  gl 
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N(Prop  ♦  Rudder) 


.  g2  ?  I/Js 
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and 
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n  *  n|  ♦  (n  -  n^ )( 1  - 

e'At/T)  |  At  =  t  -  t, 

8  *  6^  ♦  J(t  -t^  )  . 

-6l«  8  <  8l 

*.  ■  f  i*  "»i 

B0  =  f  L'  HYI 

C0  *  f  >■*  HNJ 

*i  •  f  lm«; 

*1  -  5  LHVi 

c,  *  t  L*  hm: 

*2  ■  I  i* 

B2  -  |  l-  hy; 

C2  *  J  L-  HN; 

*)  ■  *. 

Bj  .  f  L»H(Yrv;  *  y;> 

CJ  '  1 

\  ■  P<**  «! 

B.  =  r  Is  HY  ' 

4  4  rrv 

c4  ■  l 

As  «=  Pd»  a2 

*5  '  T3  l'"»rrf 
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•«  *  h  l,('wl  *  >v;> 

cs  *  T3 

A7  =  f*P  a4 

0,  =  \  1.HY, 

C,  -  |  I’M, 

Ag  =  Pd-  a5 

*0  *  I  “"2 

C,  .  {  l-HNj 

y  1 lH  *j 
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A10S  ^  *7 

»10=  P^ho 

C10*  P^  co 

A,,*  Pd»  a8 
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A12*  -9 
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*13*  P**  *10 

B,j=  Pd-bj 
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Ship  motlona*  poaltlon  and  orientation  relative  to  earth-fixed  axes  are  provided  by 

u  *  u  coat-v  alnf  +  U  coa* 

°  o 

v0  *  u  alnf  +  v  coat  *UC  •*ntc 
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r  e  r 
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♦  X  (o) 


y0 

♦ 


£  v0dt  <  y0(o) 

J*  rdt  *  r(o) 
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EXPERIMENTAL  PROGRAM 

Experiments  were  made  to  evaluate  coefficients  of  the  motion  equations. 
The  program  was  divided  into  three  parts,  according  to  the  facilities  used. 
Since  the  number  and  ranges  of  variables  are  extensive,  the  test  programs 
were  limited  to  conditions  necessary  for  application  of  the  equations  to 
particular  maneuvers. 

PROGRAM  I:  Straight-Course  Experiments 

A  14-ft  captive  model  was  towed  at  various  speeds.  Propeller  speed 
was  varied  among  positive  and  negative  values,  and  rudder  deflection  was 
set  at  intervals  between  plus  and  minus  30  degrees.  A  Series  60  model  was 
used  ,  particulars  of  which  are  given  in  Table  1. 

The  propeller  was  driven  by  a  frequency-controlled  synchronous  motor, 
inteqral  with  a  propeller  thrust  and  torque  dynamometer.  The  power  source 
and  speed  controller  were  located  ashore  and  connected  electrically  to  the 
model  through  overhead  cables  via  the  towing  carriage. 

Rudder  deflections  were  set  manually,  using  a  tiller  and  protractor. 

The  three  force  components  measured  were  the  longitudinal  component 
X  (in  line  with  the  model's  longitudinal  centerline)  and  the  two  lateral 
components  and  Yg  (perpendicular  to  the  model's  centerline,  forward 

and  aft).  The  net  side  force  is  Y  =  Y.  +  Y_  .  The  net  yaw  moment  (act- 
i ng  about  a  vertical  axis  through  the  model  CG )  is  determined  by 
N  =  ^x?/^+^bx^b  *  where  F  1S  longitudinal  separation  between  the 
model  center- of-mass  and  the  force-measuring  balances  (^^>0  •  <  0)  . 

The  force  balances  are  of  the  differential  transformer  type,  capable 
of  measuring  two  orthogonal  force  components  and  a  moment  component  about 
a  central  axis.  Balance  locations  and  constraints  are  shown  in  Sketch  III 
on  the  following  page.  The  model  was  constrained  against  surge,  sway» and 
yaw  mot  ions, but  was  permitted  freedom  to  heave,  pitch  and  roll. 
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SKETCH  III 


SECTION  B-B 


BE] 

SECTION  A-A 


Tests  in  this  program  were  conducted  on  straight  path  ( r  *  =0)  with  no 
drift  angle  (v=0).  Model  speed,  u  ,  was  varied  between  0  and  5*7  feet 
per  second.  Propeller  speed  was  varied  between  ±7*65  revolutions  per  sec¬ 
ond.  No  data  are  reported  for  propeller  speeds  less  than  4.4  rps*  for  reasons 
of  scale  effect  on  forces.  Rudder  angles  were  limited  to  ±30  degrees* be¬ 
cause  of  the  effect  of  scale  difference  on  the  point  at  which  rudder  stall 
occurs  (near  30  degrees). 

Turbulence  stimulation  (of  flow  about  the  hull)  was  aided  by  a  0.03" 
inch-diameter  trip  wire  located  8.4  inches  aft  of  the  forward  perpendicular. 
It  is  unlikely  that  turbulent  flow  was  achieved  at  most  of  the  hull  speeds 
of  this  test.  However,  hull  resistance  at  low  speeds  is  not  significant 
compared  with  thrust  forces.  Furthermore,  hull  resistance  is  deducted  from 
the  total  measured  X-force  in  arriving  at  thrust  coefficient  information 
for  use  in  ship-motion  computations;  thus,  the  error  introduced  is  subse¬ 
quently  subtracted. 

Hydrodynamic  force  and  moment  data  are  listed  for  various  conditions 
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of  ship  speed,  propeller  speed,  and  rudder  angle  (u,  n,  6), in  Table  II. 
These  data  are  dimensional  for  the  l4-foot-1ong  model. 


The  data  are  reduced  to  coeff 
siontess  coefficients  are  defined  in 
forces  or  the  section  on  huM  forces 
by 

S  =  Pn  -  u 

o  pn 

Coefficient  data  are  plotted  in  Figu 


cient  form  in  Table  III.  The  dimen- 
the  section  on  propeller  and  rudder 
Apparent  slip  ratio,  Sq  ,  i s  given 

es  2,  3  and  4. 


PROGRAM  II:  Large  Drift-Angle  Experiments 

A  5-ft  captive  model  was  towed  at  various  drift  angles  through  a 
total  range  of  360  degrees.  A  Series  60  model  was  used,  particulars  of 
which  are  given  in  Table  I.  To  avoid  wall  effects,  the  tests  were  con¬ 
ducted  in  Davidson  Laboratory  Tank  No.  2,  with  the  rotating  arm  at  its 
longest  practical  radius.  Side  force  and  yaw  moment  were  measured  through 
the  entire  range  of  drift  angles.  Since  coupling  effects  of  propeller  and 
rudder  should  be  small  relative  to  hull  forces,  they  were  neglected  in  the 
test.  The  propeller  and  rudder  were  both  fixed. 

The  model  was  attached  to  the  rotating  arm  in  a  manner  similar  to 
that  used  in  conventional  rotating-arm  experiments  of  surface-ship  models.^ 
The  principal  di fference  was  that  the  flexure  plate  mounting  on  the  model 
was  made  rotatable  in  steps  of  60  degrees.  This  extended  the  drift-angle 
capability  of  the  balance  'jeyond  the  usual  ±30-degree  limits, to  any  desired 
angle.  The  normal  arrangement,  with  fixed  mounting,  constrains  the  model 
against  yaw,  surge,  sway, and  roll  mot  ions. but  permits  freedom  to  pitch  and 
heave.  The  modification  changes  the  nature  of  the  roll  and  pitch  con¬ 
straints;  but  negligible  tendency  to  pitch  or  roll  was  detected  at  the  low 
model  speeds  tested. 

Model  speed  was  selected  according  to  the  drift  angle  of  the  run. 

For  drift  angles  near  0  or  1 80  degrees,  speeds  between  1.5  and  2.5  feet 
per  second  were  used.  At  large  drift  angles,  the  speed  0.85  feet  per 
second  was  used,  v.  o  avoid  unrealistically  high  Froude  numbers  in  condi¬ 
tions  that  would  exist  only  at  very  low  ship  speeds.  Speed  effect  at 
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90  degree  drift  angle  received  special  attent  ion. The  dimensionless  side  force 

Y 

coefficient  Y*  =  -  is  plotted  against  speed  for  this  condition, in 

|  LHU3 

Figure  5* 

Measurements  comprised  X  and  Y  components  of  hydrodynamic  force 
and  N  component  of  hydrodynamic  moment.  Oimensionless  coefficients  N'  , 
X'  and  Y'  are  formed,  and  plotted  in  Figures  6,  7  and  8,  respectively. 
The  rotating-arm  facility  imposed  a  small  path  curvature  in  all  tests  of 
this  p  ogram.  This  effect  was  suppressed  by  averaging  data  taken  at  posi¬ 
tive  and  negative  drift  angles.  Averaged  data  are  represented  in  the 
figures.  In  view  of  the  symmetry  of  the  averaged  data,  drift-angle  argu¬ 
ments  in  the  figures  are  shown  only  from  0  to  1 80  degrees.  Ideally,  the 
drift  angle  experiment  would  be  conducted  on  straight  cour  se.snd  at  many  com¬ 
binations  of  path  curvature  and  drift  angle. 

PROGRAM  III:  Pure-Yaw-Rotat ion  Experiment 

A  5-ft  model  was  moved  in  pure  yaw  rotation  by  applying  a  known 
yaw  moment  via  falling  weights  and  pulleys.  The  resulting  steady  yaw 
rate  was  measured.  A  Series  60  model  was  used;  particulars  of  which  are 
given  in  Table  I. 

A  vertical  mast  was  erected  in  the  model  coincident  with  the  model's 
negative  z-axis  (through  the  CG).  The  mast  was  a  drum  having  two  fine 
cords  attached,  which,  when  wound  around  the  drum,  applied  a  pure  yaw 
moment.  The  cords  were  led  over  ball-bearing  pulleys  attached  to  the 
outer  sides  of  the  walls  of  Tank  No.  3-  Weight  holders  were  hung  from  the 
cords, and  the  system  was  arranged  to  permit  the  weights  to  hang  freely  and 
travel  downward  as  the  model  commenced  to  rotate. 

In  testing,  the  lines  were  wound  about  the  drum,  elevating  the 
weights.  When  the  model  had  been  positioned  directly  between  the  pulleys, 
the  model  was  released.  After  an  acceleration  interval,  the  time  of  one- 
half  revolution  in  yaw  was  measured.  Different  weights  permitted  a  de¬ 
termination  of  the  effect  of  yaw  rate  on  the  dimensionless  moment  co¬ 
efficient  ,  — — -  .  Results  are  plotted  in  Figure  9.  These  show 

f  L4H  ra 

regligible  change  of  coefficient  value  over  the  r  range  of  the  test. 
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COMPUTATION  OF  SIMPLE  DOCKING  MANEUVER 


General 


The  equations  of  ship  motion  (16  through  20)  were  programried  for  nu¬ 
merical  integration  by  digital  computer.  The  computer  source  program 
written  in  Kingston  Fortran  language,  is  listed  in  Appendix  B.  The  pro¬ 
gram  applies  to  calculation  of  various  turning,  stopping,  and  accelerating 
maneuvers  of  a  single-screw  ship. 

An  example  application  of  the  mathematical  model  is  made  in  the  cal¬ 
culation  of  a  simple  docking  maneuver;  namely,  the  "one-bell"  landing. 

This  maneuver  involves  a  single  engine  reversal  to  halt  a  ship  alongside 
its  berth.  The  particular  purpose  is  to  demonstrate  sensitivity  of  ship 
response  to  operating  and  design  parameters.  Conditions  of  the  maneuver 
are  idealized  in  the  sense  that  no  tugboat  or  MPD  forces  are  applied,  no 
w’’nd  acts,  and  deep-water  hydrodynamic  coefficients  are  used.  The  one-bell 
landing  is  a  legend  for  large  ships,  but  as  a  limit  maneuver  it  provides 
insight  into  the  docking  problem. 


Characteristics  of  the  study  ship  are  given  in  Table  I.  Hull  dimen¬ 
sions  approximate  the  MARINER'S,  but  are  based  on  a  Series  60  form.  Co¬ 
efficient  values  are  listed  in  Table  IV  and  values  of  other  parameters  in 
Table  V.  Conventional  Y  and  N  hydrodynamic  derivatives  are  obtained 
from  Reference  13  (Model  60). 


To  extend  the  results  of  computations  to  ships  of  other  sizes,  the 
following  scaling  laws  may  be  used:  Subscript  A  refers  to  the  study 
ship,  and  B  to  a  geometrically  similar  ship  of  different  size.  The  fac¬ 
tor  X  is  the  scale  r,  vio  based  on  length;  i.e.,  X  =  L^/L^ 


x_  =  x,  x  X 
C5  H 


examples  of  scaling  are  shown. 
0 i s  Lance : 

Angle : 

Time: 

Force : 

Moment : 


’B 


“  *A  X  ' 


*B  =  *A  * 


x'/2 


B 

N„  =  N 


XA  x  X" 
A 


B 


x  X* 


Linear  acceleration1: 
Angular  acceleration: 
Speed 

(ship  and  current): 
Angular  speed: 


A  * 

Typi cal 

• 

UB 

=  uA  x  1 

A 

• 

r  „ 

ii 

•u 

X 

y 

i 

B 

A 

u_ 

=  u  x  X^/2 

B 

A 

n 

=  nx  X->/2 
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The  correlation  of  hydrodynamic  terms  between  models  of  different 
size  (and  especially  between  model  and  ship)  is  a  major  problem  in  ship 
hydrodynamics.  For  this  reason  a  14-ft-long  model  was  used  to  measure 
forces  depending  on  rudder  actions,  and  friction  scale  effect  was  account¬ 
ed  for  in  scaling  thrust  forces.  Still*  the  main  interest  of  the  work  is 
to  extend  ship  motion  analysis  to  permit  evaluation  of  parameters  affect¬ 
ing  maneuvering  performance  at  low  speeds,  and  not  the  precise  determina¬ 
tion  of  ship  response. 

The  basic  conditions  of  the  example  maneuver  are  these: 

The  ship  approaches  a  pier  at  a  constant  velocity  of  4.3  knots,  relative 
to  the  water.  This  corresponds  to  a  propeller  speed  of  20  rpm.  Turning 
rate  and  drift  angle  are  initially  zero.  The  "full  astern"  command  is 
executed,  and  subsequent  ship  motions  and  trajectory  are  computed. 

The  calculation  is  repeated  for  several  variations  of  each  of  the 
following  parameters: 

1.  Maximum  reversed  propeller  speed,  n* 

2.  Propeller-response  time  lag*  i.e..  time  constant  T 

— * 

3.  Velocity  of  uniform  water  current.  Uc 

With  propeller  reversed  at  low  ship  speed,  rudder  deflection  was  found  to 
have  negligible  effect  on  ship  motion  and  trajectory;  hence,  it  is  not  in¬ 
cluded  as  a  parameter  in  these  computations. 

Oistance  and  time  required  to  stop  are  used  as  performance  criteria 
for  studying  sensitivity  to  parameters  I  and  2.  For  parameter  3.  the 
lateral  trajectory  characteristics,  position  and  velocity,  are  used. 

The  computation  proceeds  from  t  »  when  the  command  n*  to  propeller 
is  given.  Formally,  the  initial  conditions  of  motion  are  u  =  u(0)  , 
v  =  r  =  u  =  v=f  =  0.  Initial  heading,  0)  ,  is  arbitrary  unless  water  current 
is  specified.  Coordinates  xq(0)  and  yQ(0)  are  also  arbitrary.  Typical 
computer  output  is  shown  in  Appendix  C. 

Resul ts 

The  effects  of  maximum  reverse  propeller  speed  are  shown  in  Figure  10. 
Clearly,  variations  of  n*  about  the  base  value  of  -45  rpm  cause  signified 
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variations  in  head-reach  and  time-to-stop.  For  the  study  ship,  -45  rpm 
corresponds  to  normal  "full  astern."  Thi s  is  the  command  most  likely  to  be 
used  for  killing  headway  of  a  fully  loaded  vessel.  The  response  to  emer¬ 
gency  ful l-astern  (taken  as  -60  rpm)  is  significantly  quicker,  and  responses 
to  lower  values  such  as  half-astern  (-30  rpm)  are  unacceptably  sluggish 
for  docking.  The  curves  representing  time-and  di stance-to-stop  are  simi  lar 
in  character  under  these  conditions.  This  is  in  accord  with  the  simplified 
relations  between  constant  retarding  force.  X  ,  and  the  time  and  distance 
necessary  to  destroy  kinetic  energy.  Energy  and  motion  equations  yield 


(m-X.)uos 

5  - - Tx— 


and 


t  =  - 


(m-X. )u 
u  o 


At  low  ship  speeds,  X  is  mainly  determined  by  n2  {  hence,  both  s  and 
t  are  proportional  to  l/n2  . 


Figure  II  indicates  how  response  time  of  the  engine  affects  distance 
and  time-to-stop.  Since  the  time  constant  T  is  a  measure  of  delay  in 
achieving  n*  ,  it  dictates  delay  in  reaching  peak  deceleration.  This  af¬ 
fects  distance  covered  in  each  succeeding  increment  of  time  after  t  • 
Time-to-stop  is  not  correspondingly  affected,  however. 


The  validity  of  the  longitudinal  ship-motion  calculation  was  tested 

by  computing  stopping  maneuvers  of  the  ESSO  SUEZ  and  comparing  results 

f  Q 

with  full-scale  test  datareported  by  Hewins.  Chase,  and  Ruiz. 


Ship  trajectories  corresponding  to  three  variations  of  the  basic 
docking  maneuver  are  pictured  in  Figure  12.  Water  current  is  the  parameter, 
A  time  history  of  significant  computational  variables  for  the  no-current 
condition  is  shown  in  Figure  13-  Computed  yaw-angle  and  dr i ft-angl e  values 
develop  significantly  prior  to  ship  stopping.  Measurements  of  free  maneu¬ 
vering  are  unavailable  for  verification  of  these  results,  but  tendencies  ap¬ 
pear  to  be  correct,  based  on  experience  with  full-scale  ships.  The  general 
preference  of  ship-handlers  for  port-side-to-pier  landings  is  based  on 
these  effects  (for  right-hand  single-screw  ships). 

Trial  computations  were  made  to  discuss  the  importance  of  steady- 
state  hydrodynamic  terms  in  low-speed  transient  maneuvers.  These  showed 
that  elimination  of  all  steady-state  hull  termsdidnot  appreciably  affect  the 
longitudinal  motion  results  shown  in  Figure  13*  However,  effects  on  lateral 
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motions  were  quite  large,  roughly  doubling  the  values  of  yaw  angle,  drift 
angle,  and  lateral  displacement.  In  another  trial,  only  linear  hull  hy¬ 
drodynamic  terms  were  included.  Results  more  in  line  with  the  basic 
computation  were  then  obtained.  However,  it  is  cautioned  'hat  the  simple 
docking  maneuver  does  not  involve  large  drift  angle  or  (dimensionless) 
turning  rate  through  most  of  the  duration,  as  might  be  the  case  in  other 
low  speed  maneuvers. 

The  effect  of  uniform  water  current  on  a  slowly  moving  ship  is  im¬ 
portant.  In  Figure  12  the  effect  of  a  one-knot  current,  setting  east,  is 
compared  to  a  one-knot  current  setting  west  (relative  to  a  north-south 
pier).  The  right-hand  ship  is  set  heavily  westward  and  put  hard  against 
the  pier.  The  left-hand  ship  is  set  eastward,  and  finds  that  a  portside- 
to-pier  landing  is  impossible.  Therefore  a  starboard-side  landing  is 
shown.  Because  uniform  water  current  adds  a  constant  velocity  to  other 
ship  motions,  time-to-stop  is  very  important  in  water -current  cases. 

A  limiting  criterion  for  docking  is  contact  velocity  with  the  pier. 

y- 

Lateral  contact  velocities  of  -2.3  feet  per  second  (for  the  ship  set 
westward  in  Figure  12)  or  1.1  feet  per  second  (for  the  ship  set  eastward) 
could  result  in  costly  side-shell  damage.  A  realistic  non-uniform  current 
would  be  weaker  near  the  pier,  and  more  confused,  but  the  desirability  ot 
having  additional  lateral  control  forces  available  is  apparent. 


Combined  with  yaw  rotation,  the  lateral  velocity  at  the  stern  (at  station 
18,  on  20  stations)  is  -2.7  feet  per  second. 
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CONCLUSIONS 

Conclusions  of  this  study  are  based  on  experimental  results,  results 

of  computed  maneuvers,  and  mathematical  model  development. 

Conclusions  Based  on  Experiments 

1.  Straight-course  experiments  show  that  the  speed  of  an  ahead-turning 
propeller  strongly  influences  rudder  effectiveness  of  a  single-screw 
ship.  The  relation  may  he  expressed  in  terms  of  the  propeller  ad¬ 
vance  ratio  for  ahead  ship  motions.  When  propeller  rotation  is 
reversed,  rudder  forces  are  small,  erratic,  and  ineffective  for  con¬ 
trol  . 

2.  A  significant  side  force  is  generated  by  the  rotation  of  a  single 
propeller  behind  a  ship,  especially  for  reversed  propeller  rotation 
at  low  ship  speeds.  This  force  is  responsible  for  a  tendency  toward 
positive  (clockwise)  yaw  rotation  during  ship-stopping.  The  ratio 
of  measured  yaw  moment  to  side  force  locates  the  center  of  effort  at 
approximately  1 0%L  forward  of  the  rudder  post, for  reverse  propeller 
rotation,  and  approximately  at  the  position  of  the  rudder  for  ahead 
propeller  rotation. 

3*  Hydrodynamic  yaw  moment  resulting  from  pure  yaw  rotation  is  closely 
proportional  to  the  square  of  yaw  rate  for  a  5-foot  Series  60  model 
at  yaw  rates  between  0.04  and  0.07  sec"'. 

Conclusions  Based  on  Computed  Maneuvers 

1.  Computations  show  the  importance  of  maximum  reverse  propeller  speed, 
n*»  in  ship  stopping  from  low  speeds.  Since  hull  resistance  is  minor 
at  low  speeds,  the  stopping  force  is  mainly  determined  by  n“  squared. 

2.  The  time  response  of  the  engine  importantly  affects  the  speed  with 
which  astern  thrust,  hence  ship  deceleration,  is  established.  Since 
this  a f fects  the  duration  of  the  high-speed  part  of  the  maneuver,  it  has 
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much  greater  effect  on  distance  than  on  total  time  required  to  stop. 

3*  For  some  maneuvers,  time-to-stop  is  more  important  than  distance.  For 
example,  ship  translation  due  to  a  uniform  water  current  depends  di¬ 
rectly  on  total  elapsed  time,  and  not  at  all  on  distance  traveled 
through  the  water.  Therefore,  a  maneuver  which  may  be  simple  for  an 
unassisted  ship  in  still  water  could  be  impossible  in  a  modest  cur¬ 
rent,  depending  on  the  time  it  takes  to  complete. 

Conclusions  Based  on  Development  of  Mathematical  Model 

1.  Propeller  forces  are  conveniently  represented  in  motion  equations  by 

followinq  the  familiar  vs.  J  notation.  Extension  of  the  k.  , 
a  t  s  t 

concept  assists  in  the  representation  of  longitudinal  force, 
lateral  force,  and  yaw  moment,  as  convenient  functions  of  ship  speed, 
propeller  rotation,  and  rudder  angle.  Applied  to  both  ahead  and  re¬ 
verse  propeller  rotation,  the  technique  provides  insight  into  the 
changing  propeller  flow  conditions.  After  blade  stall,  a  locked  pro¬ 
peller  drag  function  may  suffice  until  stall  recovery.  Computations 
may  be  organized  on  the  basis  of  instantaneous  value. 

2.  If  the  water  current  can  be  assumed  uniform,  vector  addition  of  ship 
motion  through  the  water  and  motion  of  the  water  current  permits  a 
considerable  simplification  of  the  equations  of  motion.  On  the  other 
hand,  if  water  current  is  so  irregular  as  to  disallow  this  assumption, 
mathematical  representation  may  be  impractical. 

3.  In  general,  the  variety  of  ship  motions  possible  at  low  ship  speed  is 
boundless.  Because  of  this,  appropriate  force  and  moment  data  are 
necessary  in  computations  of  low-speed  ship  maneuvers.  Steady-state 
hydrodynamic  reactions  are  complex  functions  which  cannot  be  arbi¬ 
trarily  excluded  from  ship-motion  equations,  even  at  low  ship  speeds. 

The  mathematical  model  is  a  useful  tool  which  may  be  extended  to  in¬ 
clude  additional  forces,  such  as  wi na  and  auxiliary  propulsion  devices. 
Constraints  due  to  anchor  or  mooring  lines  may  be  added.  With  the  proper 
data  on  external  forces,  the  mathematical  model  provides  an  excellent 
means  for  analyzing  many  complex  problems  associated  with  low-speed  ship 
maneuvers . 
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TABLE  I 


MODEL  AND  SHIP  CHARACTERISTICS 


HULL 

14-FT  MODEL 

5-FT  MODEL 

SHIP 

Length,  BP,  L  .  ft 

14.00 

5.00 

528.5 

Breadth,  B  ,  ft 

1.868 

0.714 

75.5 

Draft,  H  ,  ft 

0.747 

0.267 

28.2 

Block  Coefficient,  C, 

0.60 

0.60 

0.60 

l/b 

7-50 

7.0 

7.0 

L/H 

18.75 

18.75 

18.75 

Displacement,  A 

731  1b  F.W. 

35-6  1b  F.W. 

19*300  L.T. 

L  x  H  ,  sq  ft 

10.46 

1.335 

14,900 

Lateral  Area  of  Hull  ,  sq  ft 

10.35 

1 .320 

14,750 

Mass  of  Model ,  M  ,  slugs 

22.75 

1.11 

1 ,345.000 

LCG/L  from  FP 

0.515 

0.515 

0.515 

Radius  of  Gyration  in  Air,  ft 

132.1 

RUDDER 

Area  of  Rudder,  Ap ,  sq  ft 

0.167 

0.021 

238 

arAh 

0.016 

0.016 

0.016 

Aspect  Ratio  of  Rudder 

1 .90 

1  .90 

1.90 

PROPELLER 

Troost  -  4  blades 

Diameter,  d  ,  ft 

0.518 

0.198 

19.54 

Pitch,  p  ,  ft 

0.599 

0.187 

22.50 

p/d 

1.15 

0.95 

1.15 

Mean  Width  Ratio,  MWR 

0.261 

- 

0.261 

Expanded  Area  Ratio,  EAR 

0.550 

- 

0.550 

Blade  Thickness  Fraction, Btf 

0.045 

- 

0.045 

Rake,  deg 

6 

- 

6 
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TABLE  II 

STRAIGHT 

-COURSE 

FORCE  AND 

MOMENT 

DATA 

RUN  NO  . 

U 

n 

5 

X 

Y 

N 

fps 

rps 

deq 

1b 

1b 

Ib-ft 

1 03 

•  OO 

7  •  6b 

o 

4*14 

-.06 

•  50 

IUI 

•  oo 

7.65 

1  o 

4.14 

.41 

-2.V6 

vv 

•  oo 

7.65 

20 

3.  V4 

.77 

-5.56 

0V7 

•  oo 

7.60 

50 

3.64 

.  VV 

-7.06 

tUb 

•  oo 

7  •  65 

-io 

4.20 

—  .  62 

4.24 

1  u7 

•  oo 

7.65 

-20 

4.00 

-.  V6 

6.50 

iuv 

•  oo 

7.65 

-30 

3.60 

-1.17 

7.76 

71 

•  06 

7.65 

O 

3.50 

-.07 

•  45 

vo 

•  06 

7  •  6b 

1  O 

3.54 

.61 

-4.12 

V2 

•  06 

V  •  ob 

2o 

3.3  3 

1  •  06 

-7.72 

V4 

•  05 

7  •  6b 

30 

3.03 

1  .46 

-10.50 

4  3 

•  54 

4.61 

3o 

•  VI 

•  33 

—  2 .64 

72 

•  06 

7  •  6b 

-  1  o 

3.64 

-.70 

4.73 

70 

•  06 

7.70 

-20 

3.33 

-1  .25 

6.34 

04 

•  06 

7.67 

-30 

2.63 

-1  •  6  V 

10.56 

44 

•  60 

4.5V 

30 

•  OV 

.31 

-2  .VO 

64 

•  06 

4»bV 

O 

1.11 

-.06 

•  2  V 

bO 

•  OV 

4  •  6  1 

1  U 

1.05 

•  25 

-  1  .60 

52 

•  06 

4  •  5b 

20 

1.0  1 

•  44 

-3.02 

bO 

•  04 

4  •  52 

30 

.61 

•  4  V 

-3. VI 

4b 

•  OO 

4.5V 

30 

•  6  V 

.  33 

-3.17 

7V 

•  06 

4.5V 

-20 

1.0  1 

—  •  4  V 

5*61 

Ob 

•  06 

4  •  5  V 

-JO 

.6  1 

-.66 

4  .22 

1  U4 

•  OO 

-  7  •  6b 

O 

—2 .75 

-.52 

5.44 

1  02 

•  OO 

-  7  •  6b 

IO 

-2.7b 

—  .  62 

3.46 

luu 

•  oo 

—  7  »6b 

20 

-2.73 

—  .  66 

3.V6 

VO 

•  oo 

v  —  7 .65 

30 

-2.63 

-.60 

5.61 

106 

•  oo 

“7  »6b 

-IO 

-2.64 

-.27 

3.65 

tuo 

•  OU 

-7,65 

-20 

-2.63 

-.64 

3.1V 

1  to 

•  oo 

—  7 .65 

-30 

—2 .61 

-.36 

4.52 

7U 

•  06 

—  7  .  6b 

O 

-3.4  3 

—  .  6  V 

3.01 

VI 

•  06 

-7.65 

to 

—  3.33 

-1.23 

4  -62 

V3 

•  06 

-7.65 

20 

-3.23 

—  .61 

2 .27 

Vb 

•  05 

—  7  •  65 

30 

-3. 4J 

-.66 

3.V6 

77 

•  06 

—  7  •  66 

-10 

-3. 03 

-.34 

1  .VO 

63 

•  06 

-7.67 

-20 

-2.V3 

—  .61 

3.17 

OV 

•  06 

—  7  •  65 

-JO 

-3.03 

-.74 

3.41 

6V 

•  06 

-4.56 

o 

-1.21 

-.25 

1  .25 

63 

•  06 

-4.5V 

IO 

-1 .31 

—  .15 

•  66 

57 

•  06 

-4.52 

20 

-l  .31 

-.35 

1  .24 

bl 

•  05 

-4 .6  1 

30 

-1.21 

-.37 

1.17 

76 

•  06 

-4.5V 

-to 

-l.l  1 

-.35 

1  .53 

02 

•  06 

-4.62 

-20 

—  1  .  2  1 

—  .13 

•  V6 

oo 

•  06 

-4  •  bV 

-30 

-1.21 

—  .  1 6 

.65 

lb2 

4*20 

7  •  66 

O 

-.61 

-.04 

•  O  1 

154 

2  •  7  V 

4  •  V6 

0 

-.30 

-.03 

-.11 

15V 

4*22 

7.42 

O 

-1.11 

-.30 

»  JV 

16U 

4*23 

7*44 

0 

—  •  V7 

-.OV 

•  46 

165 

4.23 

7.42 

0 

-1.01 

-.16 

•  23 

151 

4.20 

7.35 

0 

-.VI 

-.03 

•  1  1 

167 

2.7V 

4.76 

0 

-.30 

-.04 

•  23 

15b 

5*65 

7.54 

o 

-3.64 

-.11 

-  .22 

161 

4*23 

4.75 

0 

-2.32 

-.  15 

-.14 

lb 

2.7V 

4*53 

o 

-.51 

-.OV 

.05 

156 

5*65 

4  .V  1 

o 

-5 . 25 

-.IV 

—  .63 

16V 

1  *44 

-4.4V 

o 

-1.11 

-.07 

.  IO 

160 

2  •  7  V 

-7.42 

0 

-2.71 

-.25 

•  33 

164 

4*23 

-7.47 

0 

—  5 .47 

.25 

-2.4 

1 66 

2.7V 

—  4  •  5  1 

o 

-2.12 

•  O  1 

-.55 

163 

4.23 

-4.4  1 

U 

-4.55 

-•OV 

IO 

150 

5*65 

-4.45 

0 

-6.6V 

-.26 

•  66 
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TABLE  III  STRAIGHT-COURSE  DATA  IN  COEFFICIENT  FORM 


RUN  NO. 

Js 

Kx 

S 

kn 

s 

o 

Y' 

N' 

103 

•  oo 

•  bOb 

-•Ol 

•  1  3 

1  .ooo 

lot 

•  uo 

•  buU 

•  Ob 

-.71 

1  *000 

vv 

•  UO 

•  4tlJ 

•  1  O 

—  1 . 33 

1  .UOO 

V  7 

•  oo 

•  4b  3 

•  13 

-1.70 

1  .OOO 

lob 

•  oo 

•  b  1  6 

-  •  ub 

1  .  O  1 

1  .OOO 

107 

•  oo 

•  40  1 

-•  13 

1  .  b4 

1  .ooo 

lov 

•  oo 

•  466 

-  •  1  4 

1 . 04 

1  .UOO 

71 

•  ^<^ 

•  44  3 

-  »o  I 

.  1  1 

•  b  1  3 

—  •  0  1 

•  uo 

VO 

•<;<: 

-•47 

.07 

-  .Ob 

•  b  1 3 

•  ob 

-•04 

V3 

•  d^ 

•  4m. 

•  1  3 

-1.63 

•  b  1  3 

•  1  4 

-.07 

S»4 

•  3  l 

•  JU'j 

•  lb 

-3  •  44 

•  bib 

.30 

-•  1  U 

43 

.3  3 

•  3  10 

•  1  1 

-1.73 

•  004 

•  1  1 

-.06 

73 

tci. 

•  4(iU 

-  »oo 

1.13 

•  013 

-.UO 

•  Ob 

7b 

•  <;<: 

•  4  1  b 

-.1*3 

1  •  7b 

•  0  13 

-.17 

•  Ob 

04 

•  33 

•  3b  O 

*•1  1 

3.40 

•  0  13 

-•£3 

•  1  o 

44 

•  3V 

•  3c  3 

.  1  1 

-  1  •  O  1 

•  7b3 

•  07 

-•04 

64 

•  36 

•  4  1  6 

-  lO  J 

•  1  O 

•  00  7 

-•Ol 

.  oo 

bb 

•  37 

•  3‘V  3 

•  ob 

-l.lb 

..6  7  7 

.03 

-.03 

b3 

•  36 

•  300 

•  1  b 

—  3 .03 

.  606 

•  06 

-.03 

bO 

•  36 

•  3<-  3 

«  1  7 

-3 .06 

•  600 

•  07 

-.04 

4b 

•  37 

•  4.73 

.  1  1 

-3.00 

.  603 

.04 

-•0  3 

70 

•  36 

•  303 

-•If 

3  •  3b 

•  607 

-.07 

.03 

bb 

•  36 

•  3  1  3 

—  •33 

4.70 

•  60  7 

-•OO 

.04 

104 

•  OO 

—  •  33  7 

-  •  06 

.03 

1  .OOO 

103 

•  OO 

-•337 

-•Ob 

.03 

1  .UOO 

IOO 

•  OO 

-•33b 

-•Ob 

.04 

1  »UUU 

08 

•  OO 

-•347 

-  .07 

•  VO 

1  .UOO 

106 

•  OO 

—  •  3<-b 

-•03 

•  8  7 

1  .uou 

108 

•  OO 

—  •  333 

-•UO 

•  76 

1  »uoo 

no 

•  OO 

—  •330 

-  »Ub 

1  >03 

1  *000 

70 

-•<;«; 

-  .400 

-•UO 

•  71 

1.16b 

-•OO 

•  03 

VI 

-  •  306 

-•  lb 

1.10 

1  •  Ibb 

-•lb 

•  04 

V3 

-••it. 

-•  JUJ 

-.07 

•  b4 

I  •  1  bb 

-•Ob 

•  03 

V5 

-•til 

-  *400 

-.11 

.04 

1  •  lbb 

-.13 

•  04 

77 

-••iti 

—  •  3bO 

-  .04 

•  4b 

1.107 

-•Ob 

•  03 

83 

-•titi 

-•34b 

-  .07 

•  7b 

1  •  IU7 

-•Ob 

•  03 

bV 

-•33 

—  •  3bb 

-  .OO 

•  bl 

1  •  lbb 

-•  lo 

•  03 

6V 

—  •  36 

-•300 

-  »ov 

•  03 

1  •  Jib 

-•03 

•  Ol 

63 

—  •  36 

-  •  4  1  O 

-  •  Ub 

•  bb 

1  .31  3 

-•03 

•  O  1 

b7 

-•37 

—  •  4  t  J 

-.13 

.04 

1  *316 

-•Ob 

•  Ol 

bl 

—  •  36 

—  •373 

-.13 

•  7b 

I  •  306 

-•Ob 

•  Ol 

7b 

—  •  36 

-•341 

-.13 

1  .OO 

1  .313 

-•Ob 

•  o  1 

b.' 

—  •  36 

-.3  7  1 

-  .04 

.63 

1*311 

-•03 

•  o  1 

bb 

—  •  36 

-•376 

-  .06 

•  bb 

1  .313 

-•03 

•  ol 

ib3 

1  *06 

•  1  03 

•  UO 

•  OO 

•  Obb 

•  UO 

•  oo 

1  b4 

1  »ob 

•  300 

-.Ol 

—  .  06 

•  063 

•  OO 

•  oo 

IbO 

1*10 

•  1  bb 

-  »U4 

•  IO 

•  UbU 

•  oo 

•  oo 

160 

1*10 

•  1  7b 

-•Ol 

.  13 

•  UbU 

•  Oo 

•  oo 

16b 

1*10 

•  1  71 

-.03 

•  06 

•  04  b 

•  oo 

•  oo 

lb  1 

1*10 

•  1  bb 

•  UO 

.03 

•  046 

•  UO 

•  oo 

167 

1*13 

•  33b 

-•Ol 

•  l  4 

•  03  1 

•  oo 

•  oo 

1  bb 

1  «4b 

•  007 

-.01 

-.  7b 

—  •  3b3 
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1 6  1 

1  •  7c 
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—  •  Ob 

-•  JO 

-•406 

•  OU 
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IbJ 

1*10 

•  1  oo 
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1  b6 

<:•<:<: 

-  •  3  1  4 

-•06 

-.36 

-•030 

•  oo 

•  oo 

160 

—  •  6  ti 

-•300 

-iUj 

.0  7 

1  •  b34 

•  Oo 

•  oo 

1 6b 

-•  73 

—  •331 

-•UJ 

•  Ob 

1  ibi.lt 

•  oo 

•  oo 

I  64 

-1  •  OO 

-•406 

•  0  3 

-  .63 

1  •  04  b 

.  oo 

•  OO' 

166 

-  1  ICO 

-•300 

•  UO 

-•3b 

3  •  034 

•  oo 

•  OU 

163 

- 1  «bb 

-  *034 

-  .03 

-.07 

3.603 

•  oo 

•  uo 

1  bb 

-3»4b 
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-.OO 
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TABLE  IV 


COEFFICIENT  VALUES  FOR  SHIP  IN  SIMULATED  MANEUVERS 


XI 

u 

- . 0068 

*3 

-  .0062 

X' 

o 

-.0075 

b0 

-  .010 

X  • 
vr 

.168 

b. 

-  .0326 

ao 

.494 

bp 

.294 

ai 

-.231 

b3 

•397 

a2 

O 

• 

1 

b4 

.12 

a3 

-.141 

bs 

.0127 

a* 

-.322 

bg 

-  .080 

as 

-.059 

N 1 
r 

-  .0107 

a0 

— 

N' 

V 

-  .095 

a7 

.300 

N' 

r 

-  .070 

3a 

.600 

1/2(N  '  +N '  ) 

'  '  rvv  r 

-1.43 

a9 

U"\ 

CM 

• 

1 

1/2  N  • 
rrv 

-  .127 

ai  o 

.500 

1/6  N  » 
rrr 

-  .304 

an 

1.50 

l/MNvv^) 

-  .061 

aia 

1.00 

Ni 

-  .095 

Y! 

V 

-.168 

-  .169 

-.305 

n3 

-  .525 

Y' 

r 

.089 

co 

.15 

1/2(Y  ' +Y '  ) 

rvv  r 

-.614 

ci 

-3.85 

1/2  Y  ' 
rrv 

-1.15 

C2 

-4.41 

1/6  Y  • 
rrr 

-3.20 

C3 

-2.62 

l/6( Y  ' +3Y ' ) 

.101 

C. 

-  .0061 

vvv  V 

-.305 

c5 

•  79 

-.65 
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TABLE  V 

FIXED  PARAMETERS  IN  EXAMPLE  COMPUTATIONS 
r ^  -  1 .UOO 

i\  =  -262  (15°) 

=  .845  fps 

gx  =  .376 

g?  =  -  .158 

q  3  =  -1.00 

At  =  1 . 0  sec(iteration  time) 
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FIGURE  I.  REFERENCE  FRAMES 


FIGURE  3.  VARIATION  OF  LONGITUDINAL  FORCE  COEFFICIENT  kx,  WITH  RUDDER 

ANGLE  AND  ADVANCE  RATIO  OF  PROPELLER  IN  BEHIND- SHIP 
CONDITION,  Js*u/nd  (u>o,n<o) 


P/Z  LHU 


FIGURE  5.  EFFECT  OF  MODEL  SPEED  ON  SIDE-FORCE  COEFFICIENT,  Y* 


P/2 L  HV 


-1.5 


-2.5 


-3.0 


FIGURE  10. 


n\  RPM 
-2.0 

n*/n  NORMAL 


EFFECT  OF  REVERSED  PROPELLER-  SPEED  ON  DISTANCE 
AND  TIME  REQUIRED  TO  STOP  SHIP  FROM  4.3  KNOTS 
(COMPUTATION  OF  SIMPLE  OOCKIN6  MANEUVER) 


Oflirr  ANGLE  $  EFFECTIVE  THRUST.  XT  ,  LBS 

RESISTANCE  ( INCL.  CENT.  FORCE  )XU  .  LBS 


APPENDIX  A 


Drag  of  Locked  Propeller 

In  the  operating  region  between  4  and  6  of  Sketch  I,  page12>  the 
propeller  force  is  essentially  a  drag  force  associated  with  propeller 
blades  continuously  changing  orientation  and  angle  of  attack  with  respect 
to  local  flow.  In  the  absence  of  detailed  experimental  data*  a  simple 
drag  coefficient  for  a  locked  propeller  can  suffice.  This  may  be  estimated 
on  the  basis  of  developed  blade  area.  DAR.  ship  speed,  u  .  and  mean  wake 
fraction  in  way  of  the  propeller,  w  . 

X(Pr°P*  )jraggjng  =  CD  5  ^  UP3 

=  C„  j  x  DAR  [u(l-w)]S 

-  pd8ae  u8 

where 

aB  =  Cp  x  DAR  x  j  x  (1-w)a 

Estimations  of  Cp  for  ship  propellers  are  discussed  in  Reference  19» 
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APPENDIX  B 

COMPUTER  PROGRAM  FOR  SIMPLE  ENGINE- AND-RUODER  MANEUVERS 

KINGSTON  FORTRAN  2.  1620.  STEVENS 
DIMENSION  C< 12.71 

I  READ  100. < ICC  I . J» « J=1 .7) . 1  =  1 ♦ 121 
PUNCH  100 

PUNCH  IOI.<(C(I.J).J= 1.7). 1=1.12) 

PUNCH  102 

CONLI  =C.<  I  .3)  »C<  I  .4)*C<  I  .51/2.0 
CONL 2  =  CONI..  I  * C  <  1.4) 


CONL3=CON|.2*C<  1  .4) 

CONL  4  =  C  ONL  0  *  C (  1.4) 

CONI  ?  =  CONL4«C< 1.4) 
COND I =C ( I  . 3)«C<  I  .6) 
COND2  =COND I «  C (  1 «  6 ) 
C0ND3=C0ND2«C ( I .6) 
C0ND4 =C0N03* C ( 1.6) 
COND5  =  COND4»C.  <  1.6) 

AO  =C (1.1 ) -CONL2  »C  <  3 . 1  ) 
A  1 =CONL I *C( 3.2) 
A2=C(I.' ) +C0NL2*C (3.3) 


A  3  =  CONP4  *  C  (  3  *  4 ) 
A4sC0ND3*C< 3.5) 
AKiaC0Nn2*CI  3.6) 
A6=C0ND4»C<3.7) 

A  7  e  C  OND  3  «  C I  4  « 1  ) 
A0=r0ND2*C( 4 .2 ) 
A9  =  C.0NL1*C(4.3> 

A  1  0  =  C0NL)4»C<  4 .4  i 
A) 1 =C0ND3*C<4.5) 
A  1 2  =  C0ND2*C (4.6) 


A  1 3  =  COND 1 *C<4.7) 
A  1  4  =(.OND4  *C  <  9  •  6  ) 


A  1  5 r- C0ND3* C  (9.7) 

HO  =  C<  1  .  1  )-C0NL2»C<5. 1  ) 
0  1  =C.ONL  1  *C  <  *>.2  ) 
B2=CONL2*C<5.3)-C<  1.1) 


D3=C0NL2*C< 5.4 ) 
B4=CONL3«C<  5.5) 
H5=CONLl *  C ( 5 .6 ) 

B6  =  C0NL4  *  C  <  5  «  7  ) 
H7=C0NL 1 *C<6. 1 ) 
t'B  =  CONL  1  *C(6.2) 

F'9  =  C0NL2*  C  <  6 . 3  )  -C  <  1.1) 
P10  =  C.0N04*C(6.4  ) 

Hi  I  =C0ND2«C.<6.5  1 


H 1 2=00NP4  *  C  <  6 . 6 ) 
B13=COND3*C(6,7> 


H 1 4  =  COND?  *C  <  7  «  1  ) 

B 1 5  =  CONL 1  *C ( 7. 2 ) 

016* C OND4  * C ( 7 . 3 ) 

C0=C( 1 .2)-CONL4»CC7.4) 
C 1 rC0NL2*C<  7.5 ) 

C  2  =  C  ONL  3  *  C ( 7 . 6 ) 

C3  =  C0NL3*C«  7.7) 
C4=CONL4*C ( P . 1 ) 
C^=CONL2*C(0.2) 


C  6  =  C  ONL  5  *  C  (  B  «  3  ) 
C7  =  C0Nl  ?*C  (  P.4  ) 


B! 


- 


J 


cn*coNL?*c<n»o» 

C9«C0NL4*C<fl«6) 

CI0»CONO5»C<B.7) 

Cl  I  «C0ND5*C<9 . 1  ) 

CI2*r.0Nr>4#C<9«2» 

CI3»CdN03«C<9.3> 

C I  4  sC0NL2*C  19.4) 

C15»C0ND5*CI9.5) 

YMPO=C< 12*1) 

FNMPn  =  YMPr>«C<  12.2)  *C  (  1*41 
J=1 

T  =  C<  1  I .5) 

u=cu r«  i ) 

V=C( 10.2) 

UG*C  <  10.1) 

VG*C1 10.2) 

R=C( 10.3) 

PPS0*C< 10.5) 

DFL=C< 10.6) 

DELO=C< 10,7) 

OELDOT  =  C (  11,1) 

X  =  C (  11  .2) 

Y=C( 11,3) 

S=C(?,7| 

PS 1=0(11.4) 

7  W=SQRTE <U«U+V»V > 

RPR  I  =  P«C<  1 .4 )/W 
SRPR1  =  ALiSE  (RPR  I  ) 

HE  T  A  =  -ATANE { V/U ) 

SHE  T  A=AHSE  I HLTA) 

RPS-r.  (  1  0  «  4  )  4  <  RPSD-C  (  10.4)  )•<  1  .0-EXPE  <-(  T-Ct  1  I  .  5  )  )  /C  <  1  .  7  )  )  ) 

IE  (SPPRI-CI2.2) )  50.50,51 

50  IE  ( SHE  TA-C (2,3))  52. 52.51 

52  IE  ( U“C (2,4) )  51,53.53 

53  XHULL =A I «U»W+A2»R«V 

YHULU=B1 «U»V4H2*U»R4B3*R*V**2/'U4B4*V*R*»2/U4R5«V**3/U+B6*R»«3/U 
ENHUL  =C I »U«V4C2»U*R4C3»R*V**2/U4C4*V*R**2/U4C5*V»«3/U4C6*R**3/U 
GO  TO  3 

51  XHULL  =A 1 «U*W+A2*R« V 

YHULL  =  B7*  V  * W4H8  * V*  ABSE ( V ) +H9*U*R 

ENHUL=C7«U*V+C0*V«U»AUSF< V*U)/W*«2+C9«R*ABSE (R ) 

3  RJ=PPS»C< I ,6)/U 

IE  (RJ-.376)  60.60,61 

60  IE  CRJ4.I50)  63,63,62 

63  IE  IRJ41.0)  65,65,64 

6 1  XC0NN=A3*RPS*«2+A4«U*RPS+A5*U**2+I A6*RPS««24A7#U*RPS) »DEL««2 
YC0NN  =  H10«RPS*»2+H1  1 «U**2+  <B12*RPS*«2+B13*U*RPS4B1 4*U#«2)«DEL 
ENC0N=C10»RPS**2+(C1 1 »RPS» »2+C 1 2*U*RPS+C 1 3*U**2 ) «DEL 

GO  TO  4 

62  XC0NN=AB*U*W4A9«U« *2»DEL* »2 . 

YCONN  =  P 1 5  *U*  *  2  *  DEL 
ENC0N=CI4»U*»2»DEL 

GO  TO  4 

64  XC0NN=A14*RPS*«2+AI5«U»PPS 
YCONN=B I  6  »RPS*  *  2 

ENCON  =  C 1 5*RPS»  *  2 
GO  TO  4 

65  XC0NN  =  A10*RPS»*2  +  AI  I • U« RP5  +  A  1 2*U* «24  A  I 3*J* * 3/RPS 
YC0NN  =  H I  6*  RP5*  *  2 

tNC  ON  =  C15»RPS*  *? 
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4  UDOT  *  (  XHULL+XCONN+C  (2*5) 1  /  AO 
VOOT  «  <  YHULL+ YCONN+ YMPD  > /BO 
ROOT*  (  ENHUL+ENCON+ENMPD l/CO 
PDI»57.28#PSl 
BFDA*57.2B*UETA 
DED=57.28*DEL 

ST  *T*0# 1 
JT*ST 
JT»JT«10 
JT  =  JT+  1 

IF((JT/J)-1)  6*5*6 

5  PUNCH  I  03*  T  *PD  I  «  U  «  X  «  Y .  BE  DA  «  RPR  I  . RPS « DED • XCONN* XHULL »UDOT »UG»VG»S 

6  U=U+UDOT*C( 2* I  » 

V=V+VOOT*C( 2. 1 ) 

R=R+RD0T*C(2« 1 ) 

PSI =PSI +R*C(2* 1 ) 

UG»C (  1 1  «  6 ) *COSF ( C ( 11.7)  )+U*COSF IPS  I  )  -V*S1NF (PS  I  1 
VG  =  C<  1  1  *6)»SINE(C(  I  1  .7)  )  +  V  #  COSF ( PS  1  )+U«SlNF (PS  I  ) 

X=X+UG»C( 2* 1  ) 

Y  =  Y  +  VG«C( 2.  1  > 

S=S+W*C(2.1) 

I F ( U  — C (2*6)  )  71.70.70 

70  IF  (  ABSF  (  DEL  '  -ABSF  <  DELD )  1  40*41.41 
40  DEL  =  DtL+C<  1  1  .  1  >«C<2. I ) 

GO  TO  2 
4  1  DEL  BDELD 
2  T=  T +C ( 2 . 1  1 
J=J+1 
'  GO  TO  7 

100  FORMAT  (55H 
1  /  (  7F  l  0  •  4  >  ) 

101  FORMAT ( 12 (7E 10.3/) ) 

102  FORMAT  (42H  T  PSI  U  X  Y  BETA  /66H1  RPR | ME 

1  RPS  DELTA  THRUST  RESIST.  UDOT  UE  VE  S) 

1  03  FORMAT  (F5.0.F7. 1 *F7.2*2Fb. 1  «F7.2/lHl.2F7.3«F6.i . F 1 0 .0 » F9 .0  * 3F7 . 3 . 
1F8.  1  ) 

71  GO  TO  1 
END 


NOTE:  Program  allows  inclusion  of  additional  constant  control  forces  in 
lateral  and  longitudinal  directions,  YMPD,  NMPO,  and  XEX.  These 
were  made  zero  in  the  present  stud/. 
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